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We consider translationally invariant states of an infinite one dimensional chain of qubits or spin-^ 
particles. We maximize the entanglement shared by nearest neighbours via a variational approach 
based on finitely correlated states. We find an upper bound of nearest neighbour concurrence equal 
to C = 0.434095 which is 0.09% away from the bound Cw = 0.434467 obtained by a completely 
different procedure. The obtained state maximizing nearest neighbour entanglement seems to ap- 
proximate the maximally entangled mixed states (MEMS). Further we investigate in detail several 
other properties of the so obtained optimal state. 
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I. INTRODUCTION 

The understanding of entanglement in a multipartite 
quantum system is a central problem of contemporary 
quantum mechanics, also spreading to statistical and 
solid state physics. A multipartite quantum system can- 
not exhibit arbitrary entanglement properties, and the 
restrictions are far from being straightforward. This fact 
has several implications on the properties of spin chains 
and spin lattices, the typical subjects of statistical and 
solid state physics. This became apparent along with 
the recent developments of density matrix renormaliza- 
tion group (DMRG) techniques (consult 1] for a recent 
review) . 

In particular, if one considers a system of two subsystems 
in a maximally entangled state, neither of the subsystems 
can, of course, be entangled with anything else. In a sys- 
tem of many quantum bits (qubits, or spin-i particles) 
the limitations on the entanglement of the qubit-pairs are 
quantified by the Coffman-Kundu-Wootters (CKW) in- 
equalities , for more than three qubits this long stand- 
ing conjecture was recently proven Q. Finding a quan- 
tum state with prescribed pairwise entanglement between 
each pair is therefore not always possible, and it is a 
rather involved task Q . 

In this paper we consider an infinite one dimensional 
chain of qubits where each qubit is entangled at least 
with its two nearest neighbours. We impose the con- 
straint of translation invariance: the state should be in- 
variant under all transformations that shift each qubit 
from its original position i to i + n for some integer n. 
Our main goal is to find the maximal possible achievable 
entanglement of the nearest neighbours, and study the 
properties of the so arising chain. This problem is inter- 
esting mainly for two reasons. First, because if the bound 
on nearest neighbour entanglement is optimal, then it 
serves as a reference point for interpreting entanglement 
values obtained for real physical systems (such as the an- 



tiferromagnetic Heisenberg chain Q or lines of ions in a 
trap which could be used for quantum computations 6]). 
Second, because the results contribute to the knowledge 
on the possible structures of distributed entanglement in 
systems of (infinitely) many subsystems. 

This question was also addressed by Wootters 0- Via 
a certain procedure he succeeded in constructing such 
translation invariant entangled chains in which the max- 
imal achievable concurrence is Cw = 0.434467 which 
corresponds to a value of entanglement of formation of 
Ef = 0.284934 ebits. This value is below Cckw = l/v^, 
the limit that the CKW inequalities would allow for, in 
the case when each quantum bit is maximally entangled 
with the rest of the system, while the bipartite entangle- 
ment is restricted to the nearest neighbours. The con- 
currence Cw is conjectured to be an absolute bound, but 
this fact is not proven. A challenging question in this 
context is if one can go beyond Cw, or even reach Cckw- 

In this paper we attack the problem in a different way, 
related to the DMRG method. This latter is found to 
be a variational method in terms of the so-called matrix 
product states (MPS) 8L equivalent to a sequence of 
entanglement swappings [3- The MPSs constitute a rep- 
resentation of the pure state of a finite number of qubits. 
For each qubit of the system an auxiliary finite dimen- 
sional Hilbert space is considered. The state is described 
by projectors acting between the auxiliary space and the 
space of the system under consideration. They are very 
suitable for approximating ground states of Hamiltonians 
in a numerically efficient way where the approximation 
lies on the dimensionality of the auxiliary system. The 
infinite chain can be studied as a limit, assuming peri- 
odic boundary conditions. This attitude is equivalent to 
a different formulation, termed finitely correlated states 
(FCS). A pure state of the whole translationally invariant 
chain is encoded into a state of an auxiliary system and 
a completely positive (CP) map. In this framework, the 
density operator for any finite subset of the system can 
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be constructed by the successive application of the CP 
map, while the auxiliary system models the rest of the 
system from the point of view of quantum correlations. 
All translational invariant states can be approximated 
in that way with an accuracy which may depend on the 
dimension of the auxiliary system. 

The entanglement distribution of such chains has been 
already studied in Refs. 0, under the restriction of 
2 dimensional auxiliary systems 0| . 
We shall use a specific subset of such FCS as an ansatz to 
maximize nearest neighbour entanglement in a transla- 
tionally invariant infinite chain of qubits. This approach 
has several advantages. First, it provides us explicitly 
with a well-defined pure state on the whole chain. Sec- 
ond, we can study a class of quantum states different to 
the one studied in Ref. 0. Third, we can also investigate 
higher correlations as next nearest neighbour entangle- 
ment and entanglement of one qubit with the rest of the 
chain. 

Thus we can study the entanglement distributed along 
the chain, however, we focus mainly on maximizing near- 
est neighbour entanglement and the properties of such an 
optimal state regarding the entanglement of next near- 
est, and further neighbours. For two dimensional aux- 
iliary systems we have analytical results describing the 
properties of entanglement of such an entangled chain 
and for higher dimensions we have performed a numer- 
ical optimization. We have found that the achievable 
nearest neighbours entanglement seems to converge fast 
to Wootters' bound Cw- 

This paper is organized as follows. Section |n] is devoted 
to the construction of the states of translationally invari- 
ant finitely correlated chains utilized in this paper. In 
Section ^Q] analytical results are presented for the max- 
imization of the nearest neighbour entanglement for the 
case of low dimensional auxiliary Hilbert spaces, to give 
an insight into the nature of the problem. In Section Hvl 
our numerical results are presented for higher dimension- 
alities, while in Section |V] the properties of the so found 
optimal states are discussed. In Section IVII our results 
are summarized, and conclusions are drawn. 



II. CONSTRUCTION OF THE TRANSLATION 
INVARIANT ENTANGLED CHAIN 



Before describing the rather mathematical construction 
in detail, we outline first the idea behind it. The trick is 
to describe the part of the chain in which one is not inter- 
ested in by a density matrix on an auxiliary Hilbert space 
TLb, a local operator of the bounded operators B{TiB) 
and a completely positive map which maps tensor prod- 
ucts of such an auxiliary local operator and the local 
operator of one qubit Ai always back into B{TL]3)- In 
this way the completely positive map ensures that one 
only archives permitted states of a qubit or more qubits 
in such an infinite chain. It is clear that the set of the 
permitted and thus realized states of a qubit or more 
qubits of the chain can increase if the dimension of the 
auxiliary system describing the "rest of the chain" is in- 
creased. Further one has to carefully choose the appropri- 
ate completely positive maps in order to study the class 
of permitted one or more qubit states one is interested 
in. 

Explicit examples for choices of the completely positive 
map for increasing dimensions of the auxiliary system 
which maximize nearest neighbour entanglement (our 
main goal) are then given in Sec. Illll but let us first 
proceed with the construction of the finitely correlated 
states in more detail. 



A. Construction of finitely correlated states (FCS) 

In the following we summarize the exact mathematical 
construction of translationally invariant finitely corre- 
lated states according to Ref. [a] . We denote by Az an in- 
finite spin-chain, the spins at sites i S Z being described 
by the algebra {A)i = M2 of 2 x 2 complex matrices, 
i.e. we describe spin-i particles or generally qubits. The 
infinite algebra Aj, arises as a suitable limit of the Zo- 
caZ tensor-product algebras yl[_„.„] := ®"=_„(-^)j- 
state Lo over A-i is specified by density matrices p[i.„] 
defining the action of lo as an expectation over local op- 
erators ^[i,„] e A[i^n]- 

^{A[i,n]) = Tr[i^„] (^p[i^„] Ai «) «> • ■ • A>) 

= Tr[i_„](p[i_„] A[i^„]j . (1) 



We consider an infinite ensemble of qubits arranged along 
a line. The first question we address is what we mean by 
the word "state" as applied to infinitely many qubits. We 
adopt the standard approach described in e.g. Ref. [l3l |. 
A state oj of the infinite chain is a functional that assigns 
to every finite set of local operators — Ai®A2 ■ ■ - ^ 

An a normalized density matrix ti;(A[]^ „]) describing the 
properties of the of n qubits. Moreover one demands that 
if one considers a subset of local operators ^[i,fc] of the 
set A^i^n] then the state u){A\^i^^) has to be obtained by 
taking the partial trace of uj{A^i n\) over the qubits not 
in A[i.„]. 



The p[i^„]'s must satisfy the compatibility conditions, i.e. 
acting locally on the n -I- 1 qubit with a unity operator 
should give the same expectation value 

w(^[i,„]) = Tr„+i (^p[i^„+i] A[i_„] (g) l„+i^ 

= Tr[i^„](^P[i^„] A[i^„]^ . (2) 

Whereas, translation-invariance requires that doing no 
operation on the first qubit should also do no change to 
the expectation value, i.e. shifting the line of qubits by 
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one (generally by an integer) 

= Tr„+i(^p[i^„+i] li (g)y4[2,„+i]) 

= Tr[i,„] |^p[i,„] v4[i_„]^ . (3) 

The class of translation-invariant finitely correlated 
states (FCS) over Az is defined by a triple {B, p, E) where 
S is a 6 X 6 matrix algebra B, pB ^ B a density matrix 
and ^ : A® B ^ B a. completely positive unital map, 
which in Kraus form reads 

¥.{A®B) ^^Vj{A®B)v] , Vj : C^^C'' C'' , (4) 
j 

I 



the r.h.s. recursively defines local density matrices 
over ^[i,n] and a total state uj on Az- Translation invari- 
ance condition Eq. 10) can be formulated as 

Trs {PB ^{1 A <»B)) = TrB{pBB) VBeB. (7) 

B. Constraints of unitality and translational 
invariance 

Concretely, we choose B — Mb {b x b complex matri- 
ces) and E like in but with just one Kraus operator 
V : «) ^ C''. (Note that this restriction does 
not decrease the generality, as the introduction of fur- 
ther Kraus operators can be avoided by increasing the 
dimensionality b.) This is such that V\ai ip) = Vi\ip), 
^^IV') = X]i=i kj) ^IlV") and with the one qubit oper- 
ator A = J2 o-ij \ai){aj\ 

E(\a,){aj\(g,B^ ^v^Bv] , BeB, (8) 

where \ai^2) G are orthonormal and ui,2 bxb matrices. 
In this notation the unitality and translation invariance 
reads for the two matrices 

viv\ + V2vl = li^ unitality (9) 

2 

vl pb Vi — pb translation invariance . (10) 

i=l 

If there exists a unique p fulfilling the previous condition, 
the resulting translation-invariant FCS are pure states 



with A & A and B £ B. Unitality means that identities 
are preserved 



E(l^ (g) le) = Ib unitality. (5) 

Let E^-'^^ (A) := E(yl(8) Ig); this defines a completely posi- 
tive map from A into B. Analogously, the recursive com- 
positions E^"^ := Eo ^id^ (g)E("~^^^ are completely pos- 
itive maps from .4[i „] into B. Setting 



(6) 

I 

over Az , namely they cannot be decomposed as mix- 
tures of other states. These pure states can be interpreted 
as ground states for appropriately constructed Hamilto- 
nians of finite range Q ■ It is also shown that this class 
of states is dense in the set of all translation invariant 
states [l^. 



The two conditions Eq. 0, Eq. H1Q|I can as well be inter- 
preted in the context of open quantum systems. I.e. is 
the state of some open quantum system where the f 's are 
the operation elements for the quantum operation. These 
operation elements satisfy the well known completeness 
relation (first condition) which leads to trace conserva- 
tion of the completely positive map, i.e. no information 
of the whole system {B plus environment) is lost. The 
second equation (condition 2) can then be interpreted 
as searching for states which are invariant under these 
interactions. 



w(A[i^„]) = Tr(p[i,„]yl[i^„]) := TrB{pB E[Ai «) E[A2 » E[ • • • (g) [An<^lB] •••]]]) 

eB 
eB 

^ ^ ^ 

eB 

^ ^ ^ 

eB 

= TrB(pBE(")(A[i,„])) , 
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C. Density matrices for a subset of qubits in the 
chain 



Let us first consider the state of one qubit in tlie chain, 
defined in Eq.© 

oj{A) = Tr^(p[i] A^Ib) := TrsipB E(A® Is)) 
= Trs (pB E (^ay \a,){a.j\ (g) Is)) 
= Trg (p_B ajj VjlBv]) 

= Tr^igiB '^^-s {"]PB'"i) \at){aj \ ® Is) 

(11) 

where we have used the cyclic property of the trace - 
Tt{XY) = Tt{YX) - and the fact that the trace opera- 
tion is invariant under the map ¥ : B ^ B dual to 
E: A^B ^ B. This is defined by 



¥{B) := J2 \aj){a.^(^vlBv, 



(12) 



Hence comparing the last term with the first one we find 
that the density matrix of a qubit in the chain is given 
by 



P[i] = Pi=Yl Tr(wj PB v^) 

Tt(vIpbVi) Tr(vlpBV2) 
Tr{vlpBVi) Tr(4pBW2) 



(13) 



Now let us investigate the nearest neighbour state, i.e. 
Eq.© by setting n = 2, thus we have ^[1,2] = Ai (g) A2 
and 

Tr(p[i,2] Ai (g) := Ttb {pb e(Ai ® ¥{A2 ® 1b 



Again using the properties of the trace-operation, the 
action of E becomes the action of its dual map F 

Tr(p[i,2] ^18)^2) := Tva®b{v{pb) Ai®¥(A2®1b 



(15) 

This provides a state PA®B E(pb) — pbV on A® 
B: 

I „ t (v\pbVi v\pbV2 

\s){t\ ® vlpBVt= [ 
^j^l \V2PBV1 V^PBV2 

(16) 

which encodes the properties of all the correlations be- 
tween one site with the rest of the whole chain. 
The right hand side of Ea. (|15l) reads Tr^,g,B (^PA(S)B Ai 

E {A2 (8) Ig)^ , by turning id^ E into its dual, nearest- 
neighbours states arise as pi2 P[i,2] = Trg^id^ (g) 
¥{pa0b) I which reads 



P12 = E 1^)01"^^'^"'' 



ij=i 



R2ijl R2ij2 



(17) 



where Rijim = Tr(i;^^D] pB viVm)- 
In general local density matrices are constructed by 

Plhn] = MvlpBVt) , (18) 

s,t 

where |s) = \si (g S2 ■ ■ ■ s„), vt := Vt^ ■ ■ ■ Vt„. 



(14) 



III. MAXIMIZING NEAREST NEIGHBOUR 
ENTANGLEMENT 

A. General discussion for an optimal choice of v\ 

Our main goal is to maximize the nearest neighbour en- 
tanglement of state Eq. (|17|l 



/Tr(i;{t;| pB ViVi) 

Tv{vlvl PB ViVi) 

Tr(-y|i>| PB vwi) 
\Tr(w|t;| Pb viVi) 



Tt:{v\v\ Pb ^1^2) 
Tt:(v\v\ Pb W1U2) 

Tx{v\v2 PB V1V2) 

Tr(4w| PB V1V2) 



Tr(w{w} Pb V2V1) 
Tr(w^i;| Pb V2V1) 

Tl{v\v2 PB V2V1) 

Tr(44 Pb V2V1) 



T:r[v[vl Pb V2V2)\ 
Tt:{v\v\ Pb V2V2) 
Tt:{v\vI pb V2V2) 
Tt:{v\vI Pb V2V2)/ 



To do so we have to choose appropriate matrices vi , V2 
satisfying the condition in Eq. and then to derive a 
unique ps satisfying the condition in Eq. (|10|l . This then 
defines the nearest neighbour entanglement given by the 
density matrix above. 



What are appropriate choices for vi? 

Obviously the entanglement shared by two neighbouring 
qubits in a chain cannot be maximal, e.g. one of the 
four pure entangled Bell type states (If/*^) = i 
111)}, = -j={\01) ± |01)}). In this case every other 
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qubit has to be disentangled with that pair. Therefore 
the reduced state of the infinite chain we are looking for 
is not pure. On the other hand it should be far away 
from the tracial state (i.e. the complete mixture) as well. 
It seems to be plausible though not conclusive that the 
reduced state vanishes on some subspace. Considering 
the basis states |00), |01), |10) we notice that if we 
choose |01) states then because of translation invariance 
we need as well |10) states. Thus we expect our reduced 



state to be orthogonal to a separable pure state, e.g. |00). 

Translated to our nearest neighbour state Ea. H19(l we 
need that e.g. Tr(w|w| wiwi) vanishes. Since pB has 
to be strictly positive (otherwise we could reduce the di- 
mension of B) it follows that vi has to be nilpotent, i.e. 
vivi = 0. Then the nearest neighbour density matrix 
Ea. (|19() gets the form 



Pl2 











^ \ 

Tt{vIvI pb V1V2) Ti(vlvl pb V2V1) Ti{vlvl pB V2V2) 

Tt:{v\v\ Pb V1V2) Tt(v\v\ pb V2V1) Tx{y\v\ pB V2V2) 

VO Tr{v\v\ Pb V1V2) Tt{v\v\ pb V2V1) Tt{v\vI pb V2V2)) 
/O 
Tr(w| Pb wi) 

lh:{v\v\ Pb vi_V2) 
VO Tr(i;| Pb vx) 

/O 

^ A B C 

B* A C 

\ C* C* 1~2A 



\ 



Tt(vIvI Pb V2V1) Tt{vI pb V2 
Tr(w| Pb vi) Tt{vI pb V2 
Tr(w| Pb vi) 1 - 2Tr(t;| pB vi)/ 



(20) 



This form of the density matrix is similar to the choice in 
Ref. except that C is required to be equal zero there. 
It means that the state considered in the reference is also 
invariant under local rotation of one qubit around the 
X- and y-axis, a crucial assumption in the construction. 
As we drop this assumption, we can test another class of 
candidates within our framework. 

Let us now discuss the properties of a density matrix of 
the form Eq. (jlll A|l . Its eigenvalues are {0,A — 1 — 
(A - \B\) ± v/(3A+|B|-l)2+8|C|2} and A e [0, 1/2]. 
First we consider its purity as measured by Trp^2 which 
equals to 1 - 4A + &A^ + 2\B\^ + A\C\^ and the purity 
for the reduced matrix, i.e. the one-qubit state, equals 
to Trpl = l-2A + 2A^ + 2|Cp. Intuitively, we expect 
the entanglement to increase for a density matrix becom- 
ing purer, while the purity of the one-qubit state should 
decrease. However we will notice that for the two-qubit 
state the opposite is true. 

As a measure of entanglement we use the concurrence 
C, introduced by Hill and Wootters [THl IT^. which is a 
monotonic ally increasing function of the entanglement of 
formation [l7| . The concurrence of a density matrix p is 
given by C{p) = max{0, Ai — A2 — A3 — A4}, where \j are 
the square roots of the eigenvalues in decreasing order 
of the matrix where p = {uy ® ay)p*{ay ® Uy) and 
p* denotes complex conjugation in the standard basis. 
For the above form of density matrices the eigenvalues 
-^pp = {A -|- ^ — \B\, 0, 0} are independent of C and 
the concurrence is simply C(pi2) = 2\B\. Clearly, we 



have separability only for vanishing \B\. 

Another useful quantity is concurrence of assistance 0, 
ITgl which is defined as the sum of the square roots of 
the eigenvalues of pp, i.e. in our case its simply given 
by Cass(pi2) = 2/1. This quantity characterizes the max- 
imum entanglement of a selected pair of qubits available 
on average when the rest of the system is subjected to 
measurements. 

In order to maximize nearest neighbour entanglement we 
have to maximize the function 



C{pi 



2 \B\ 



Tr(4i'I Pb V2V1) 



(21) 



There is an additional symmetry of the nearest neighbour 
state, namely pi2 p*i2 does not change the amount of 
entanglement. For dimension 6 = 2 the only effect is to 
introduce irrelevant phase factors to Tr(w|t)| pB V2V1). 
Though for higher dimensions b this is not the case, we 
restrict ourselves to real generators as a natural choice to 
reduce the set of parameters in the calculations. We have 
checked complex extensions numerically in the region of 
the obtained maxima where we always have found that 
it only reduces the amount of entanglement. 

In the next section we analyze the case of auxiliary matri- 
ces of dimension 2 where we can give analytical solutions 
and analyze the generalization for higher dimension. 
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B. Analytical results 

For dimension b ~ 2 we can give analytical results and go 
through the whole calculation in order to understand the 
procedure in more detail and its generalization for higher 
dimensions. This will also strengthen our assumption for 
the generators of the completely unital map £. First we 
investigate the set of solutions pB for a nilpotent vi and 
its generalization for higher dimensions. Then we discuss 
the implementation to the entanglement to the nearest 
and next-nearest neighbours. 



1. Solutions for the auxiliary density matrix pB 

We have to solve the two conditions: unitality Eq. and 
translation invariance Eq. (|10|) . In 6 = 2 dimensions the 
only nilpotent matrices are the ladder operators where 
we choose without loss of generality one of them with a 
weight cosai: 



(Bloch) 



Vl 



cosai 





1 



V2 = 



1 

sinai 



cos 01 

— sin(/)i 



sm(; 
cos ( 



(22) 



where V2 is the most general (real) solution satisfying the 
unitality condition in Eq. Eq. 

Let us now discuss the solutions for the auxiliary pB, con- 
ditions in Eq. H1(J|) . The set of possible density matrices 
which are invariant under the above chosen interaction is 
illustrated in Fig.^ where we used the Bloch sphere rep- 
resentation. Every one-qubit state can be decomposed 
into three Pauli matrices 



1 



PB 



(Is 



, (Bloch) 



(Bloch) X 2 _ 



; (Bloch) 



< 1 



(23) 



For 



^(Bloch) I 



< 1 the state is mixed (corresponding to 



Trp^ < 1) whereas for 



ch)|2 ^ 



1 the state is pure 



(Tr — 1). This real three dimensional vector 7j;(Bioch) 
is called the Bloch vector and thus the state space of a 
qubit can be represented by a sphere, where the vectors 
with = 1 are pure and cover the surface of the 

sphere, inside the sphere we have all mixed states and the 
origin represents the totally mixed state, i.e. the tracial 
state. As we consider real generators of the interaction 
the y-component of n is zero and all possible one- 

qubit states are represented by the area of a circle. And 
because of the specific choice of the ladder operator only 
states in the upper half can occur as solutions. It turns 
out that the solution is an ellipse in this Bloch's sphere, 
i.e. the following equation holds for all ai, <j)i 



(Bloch) 



0)^ 



^^(Bloch) 



1)2 

2' 



(—V 



(24) 



2[p 



B 12 



and 



where the Bloch components are n\ 
= 2[pb]ii - 1. 

For 6 = 3 the auxiliary density matrix is described by 
a qutrit state which can be decomposed analogously to 
the qubit case into Gell-Mann matrices A""^, . . . , (con- 
sult the Appendix for their definitions) 



PB 



1 
3 

L 



1 



V3n,P'°^''^ A' 



; (Bloch) 



< 1 . 



(25) 



where the Bloch vector is a 8 dimensional real vector 
with similar properties as in the qubit case. Notice that 
the positivity of the state does not hold for all vectors 
in the 8 dimensional sphere. The Gell-Mann matrices 
satisfy the similar relations as the the Pauli matrices, i.e. 
TrA' = 0, TrA'A^ — 2 6^K If we choose an analogous 
nilpotent generator 



Vl = cos ai 




V2 = diag{l,smai,l).R (26) 



where i? is a 3 x 3 (real) orthogonal matrix which we can 
build up with three angles each representing a rotation 
in the two-dimensional subspace. It turns out that the 
solution of the Bloch vectors for pB again is an ellipse, 
more precisely a 5 axial ellipsoid if we fix the three angles 
and vary only ai. Varying one angle we obtain again an 
ellipsoid but with a different center and semi-axes. 
For the density matrix pB which maximizes nearest 
neighbour entanglement the length of the Bloch vector 
in 6 = 2 is \n^^^°'^'^^\ — It turns out that for higher 
dimensions b the length of the generalized Bloch vector 
is always around see Tabled The "purity" of the 
state of the rest of the chain ps measured by the squared 
length of the Bloch vector (Sioc/i) |2 _ seems to 

be quite constant when nearest neighbour entanglement 
is optimized. 



2. Properties of the nearest neighbour density matrix pi2 

Let us now return to the original question, i.e. to 
the function that maximizes nearest neighbour entan- 
glement. We have noticed in Sec lIII Ai that due to the 
nilpotent choice of vi the trace of pj5 gives only two non- 
vanishing eigenvalues, i.e. {A-|-|_B|,A — |_B|}. Concur- 
rence of nearest entanglement for dimension 6 = 2 is 
therefore 



C{pi2) = 2Tr{vlvlpBV2Vi) = 2 \B\ 

cos'^ ai (1+sinai) sin^ <pi cos"^ (f>i 

cos^ Qi cos^ (j)! ( — 1+sin ai ) — 2 (1+sin ai ) sin-^ i 

and concurrence of assistance 



cos-^ ai cos^ 



cos cti (1+sin ai) sin cpi 
61 (—1+sin ai) — 2 (1+sin ai) sin-^ 0i 



(27) 



(28) 
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FIG. 1: (Color online.) In the Figs.(a)-(c) the upper half of the x-plane of the Bloch sphere is plotted. The possible ps's and 
the nilpotent vi satisfying Eq. 12211 form the ellipse given in Eq. lHU. The dots in Fig.(a)/(b)/(c) represent the solution for 
the Bloch vectors for ai = 0/^/0. 477r, where (j>i varies G [— f, f] with a step size 7r/128. For ai = 0, Fig. (a), the solutions 
cover homogenously the ellipse and for ai n/2 they concentrate at the tracial state, i.e. there is an decrease of the set of 
possible solutions. For (j) = the pure spin up | ■(]■) is always a solution. The two big (green) dots are the solutions for which the 
nearest neighbour entanglement maximizes, this is when the Bloch vector |n'^'°'^'''| = l/\/2 which is also plotted (inner half 
circle). For higher dimensions b the absolute value of the generalized Bloch vector always around ~ l/\/2, see Table 

Q Hence maximal entanglement between nearest neighbours is obtained if the "rest of the chain" is quite equally weighted in 
its purity and mixedness, measured by the squared length of the Bloch vector. 



We have plotted both functions in Figure IIII B 21 One 
notices that while one can obtain for the concurrence 
of assistance all possible values, i.e. Cass(pi2) G [0;!]; 
concurrence of nearest neighbour entanglement has a 
maximum value of Cmax(pi2) — — 1 = 0.41421 for 
ai = 0.427079, 4>i = 0.571859. This gives a concurrence 
of assistance of Cass(pi2) = 0.585787. 
Further one notice that concurrences and concurrence of 
assistance are only equal for ai or 0i equal zero. 
Armed with this analytical experience we proceed to the 
numerical procedure and present the results for increas- 
ing dimensionality b of the auxiliary system. 



IV. NUMERICAL OPTIMIZATION OF 
NEAREST NEIGHBOUR ENTANGLEMENT 

First we discuss the parametrization and our different 
strategies to numerically maximize nearest neighbour en- 
tanglement, then we discuss the results of the maximum 
in different dimensions b. Then we proceed with a discus- 
sion of the properties of such a chain maximizing nearest 
neighbour entanglement. 



A. Parametrization for dimension b 

The triple {B, E, ps) defining the finitely correlated state 
is obtained by a finite number of parameters. We choose 



B = Mf), and carry out calculations for different dimen- 
sionalities b. The completely positive map E is described 
by the two matrices vi and V2 which are b x b matrices 
and have to satisfy the conditions Eq. ^ I10() . 



We choose vi to be a nilpotent operator as argued in the 
previous section, with a matrix of the form 



/ 

cos(ai) 



Vl 







cos(a2) 





(29) 



described by [5/2] real parameters, [. . .] denoting the inte- 
ger part. Though the parametrization is periodic with a 
periodicity of 2tt in each parameter, the parameter values 
are unconstrained, which is an important simplification 
in the case of numerical optimization. In order to satisfy 
the unitality condition in Eq. (Q, we set 
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FIG. 2; (Color online.) In Fig. (a) we have plotted the concurrence Ci2(pi2) of the nearest neighbour entanglement in Eq. I27II 
for the parameters ai G [0,7r] and G [~f i f]- The concurrence function is invariant under a shift of qi with 7r/2 and 
(j) —(j), thus we obtain four local maxima. 

In Fig. (b) the concurrence of assistance Cass(pi2) in Eq. 1281 is plotted. All possible values can be obtained. 

All plotted quantities are dimensionless. 



V2 



= diag(l, sin(ai), 1, sin(a2), . . .) R, 



(30) 



where R is an arbitrary b x b unitary matrix. However, according to our numerical experience for up to 6 = 6 
dimensions supports the conjecture that it is enough to consider real orthogonal matrices as R. The introduction of 
general unitary R-s did not lead to the increase of the maximal nearest neighbour concurrence. Thus we build up 
the generic R € SO(b) from rotations in two-dimensional subspaces, yielding the following (periodic, unconstrained) 
parametrization of R, with 6(6 — l)/2 parameters 20] 



R = 



( cos(0i) sin(0i) . . .\ 
— sin((/)i) cos((/)i) 








1 



V 



/I 

cos(0h) sin(0fc) 
-sin(0b) cos(0b) 



I cos(02) sin((/)2) ...X / cos(06_i) 

1 ... 

- sin(02) cos((/)2) ; 

V : : : \-sin(0fc„i) 
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sin(06_i) \ 




1 

cos(0b(fc-i)/2) sin(4(b_i)/2) 
-sin((/)f,(h_i)/2) cos(0b(fc-i)/2)y 



(31) 



Thus given the dimensionality 6, and a set of parame- 
ters a = ai... a[b/2], (/) — 4>i ■ ■ ■ «b(&-i)/2j we can readily 
evaluate vi{a) and V2{a,4)). Having these matrices at 
hand, we can calculate numerically p from the transla- 
tional invariance condition in Eq. (|10|l . This can be done 
by noticing that Eq. H10|) is linear in the matrix elements 
of p, thus we have to calculate the nuUspace of the linear 
mapping 

2 

m^Y.'^U^^-p^ (32) 

in the linear space of 6 x 6 matrices, in which all the 
vectors p are suitable for our aims. We have found that 
for all the parameter settings arising in our optimization 
procedure rank L = 6 — 1 holds within the numerical pre- 
cision, therefore the nuUspace is one-dimensional. Hence 



for a fixed b and parameters a, (j>, in addition to vi (a) 
and V2{a, (j)), we obtain a unique p(a, 4>). As a numerical 
check we verified that the solution is Hermitian positive 
semidefinite in all cases which have occurred. 

Performing the above calculations we can compute the 
nearest-neighbour density matrix. From this density ma- 
trix we can evaluate the concurrence. 

Thus for a fixed dimensionality 6, we have a function 
C(a, (/)) which we can numerically evaluate. This is the 
subject of an unconstrained numerical maximization in 
terms of its parameters. Unfortunately, it is not a convex 
function, thus there is no warranty to find a global maxi- 
mum numerically. In addition the function might be not 
differentiable at certain points due to the properties of 
concurrence. Therefore we chose the simulated annealing 
method, which is known to be effective for mildly non- 
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convex and non-difFerentiable function. We have used the 
routines available in the MINTOOLKIT 0| package of 
GNU Octave software . First we have searched for the 
maximum using the samin routine, a simulated anneal- 
ing code based on the implementation by Goffe We 
have set the control parameters of the routine to nt — 20, 
ns = 10 rt = 0.75, neps = 5 and eps = lO""'^" (con- 
sult the documentation .21.] of the routine for their exact 
meaning). The routine showed a normal convergence in 
each case. Then the so-obtained maxima were used as an 
initial condition for a conjugate gradient search bf gsmin, 
with numerical gradient. We have found that the func- 
tion is indeed differentiable around this maximum. The 
conjugate gradient search showed a strong convergence. 
The so obtainable final result is somewhat more accurate 
than the one obtained directly from simulated anneal- 
ing. As a result of these procedure, we have obtained 
the parameter sets a, (p for which the nearest-neighbour 
concurrence C{a,4)) has a maximum value. Though this 
procedure does not give a full warranty for finding the 
global maximum, it is very likely that the obtained max- 
ima are indeed global. 



B. Numerical results of the maximum nearest 
neighbour entanglement 

We have summarized the results of the above described 
optimization procedure in Table In case of entangled 
chains it is conjectured that the maximum value of 
nearest neighbour entanglement as measured by concur- 
rence is Cw = 0.434467. As it is apparent from the results 
m Table n in our framework we can obtain a state which 
almost reaches this upper bound. Thus the translation- 
ally invariant finitely correlated chains can approach the 
state of an entangled chain with maximal bipartite en- 
tanglement quite fast. This is our main result. The ap- 
proximation improves with the increasing dimensionality 
b of the auxiliary Hilbert-space B. 

In addition, with accidental conjugate gradient searches 
we could obtain local maxima C = 0.43406 for 6 = 8 and 
C = 0.434095 for 6 = 9. These both constitute about 
0.09% relative difference from Cw = 0.434467. 
We have to remark however, that as the maximum value 
of Cw is a conjecture, too, and we can neither fully war- 
rant the global maximum, nor check the b oo case in 
the numerical framework, we cannot exclude the possibil- 
ity to go beyond Cw ■ Nevertheless we can prove explicitly 
that the bound can indeed be obtained in the framework 
of FCS, even under several restrictions. 
Further we have checked numerically the possibility of 
using unitary instead of orthogonal matrices for R, and 
also the application of a more general, non-nilpotent vi 
by adding certain elements to its upper diagonal. We 
have found for 6 up to 6 that this does not improve the 
obtained maximal concurrence. In addition, the so aris- 
ing vi was always nilpotent, with numerically the same 
matrix elements as in Eq. H29|) , though eventually ordered 



Concurrence 




0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 3: (Color online) Here concurrence versus purity mea- 
sured by Trp^ is plotted. The points below the curve rep- 
resent all possible two-qubit density matrices. The curve it- 
self are the maximally entangled mixed states (MEMS), first 
introduced in Ref. The MEMS are defined up to lo- 

cal unitary transformations by p^EMS = (i + |)|^+)(,/,+ | + 
(i - + for q G [0,|] (green) and 

^MEMS ^ + (1 _ g)|ll)(ll| for q G [2,1] (red). 

The dotted curve below (pink) are the density matrices of 
the isotropic states, the so-called Werner states pwornor = 
l^l+p\il!'^){il!'^\. The horizontal dashed line represents 
all states having a concurrence of Cw ~ 0.434467. The dots 
are the nearest neighbour density matrices maximizing near- 
est neighbour entanglement where from left to right the di- 
mensionality b of the auxiliary system increases. One notices 
that the solutions seem to approach the bound of the real- 
izable states, in particular the MEMS states. Further note 
that the picture of entanglement measure versus mixedness 
strongly depends on the chosen measure of entanglement and 
mixedness and that different entanglement measures define 
different ordering of states, see Ref. [2,'ii2fij. 



in a different form in the matrix. This supports the as- 
sumptions we have made as well as those in Ref. 0- 

As we explicitly can calculate elements of the nearest 
neighbour state (see Table lll]> . we can ask which final 
state is approached for enlarging the dimensionality b 
of the auxiliary system. For this we plotted the ob- 
tained nearest neighbour density matrices maximizing 
entanglement in a concurrence versus purity diagram 
(see Fig. I^J. For dimension 6 = 9 the relative differ- 
ence of concurrence and purity of the nearest neigh- 
bour state maximizing entanglement and the maximally 
entangled mixed state (MEMS) Ref. 0], in our case 

(i + iCw)lV'+>(V'+l + (I - iCw)lV'-)(V'-| + 

is concerning concurrence 0.09% and concerning purity 
1.6%. 
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Dimensionality b 


2 


3 


4 


5 


6 


7 


Nearest-neighbour concurrence 


0.41421 


0.41825 


0.43200 


0.43247 


0.43336 


0.43381 


Relative difference 

from Wootters' bound (%) 


4.66 


3.73 


0.57 


0.46 


0.25 


0.15 




0.427079 


3.27378 


0.252679 


6.345324 


3.84312 


2.71122 


Q!2 






2.888910 


0.269592 


0.10177 


3.14860 


as 










3.10541 


3.29590 


01 


0.571859 


3.14062 


0.062823 


6.22996 


5.88873 


6.27750 


02 




0.56623 


5.504548 


2.351162 


6.10731 


2.50188 


03 




4.17472 


5.892460 


2.713085 


1.48352 


3.33956 


04 






0.805037 


0.047930 


4.71882 


6.25125 


05 






0.272233 


5.137121 


1.38430 


5.62825 


06 






0.741237 


0.417055 


0.79196 


3.76442 


07 








5.628356 


4.81583 


1.09039 


08 








1.759880 


2.01345 


3.43100 


09 








5.728579 


0.306965 


3.23516 


010 








1.193187 


5.68444 


2.87925 


011 










6.03621 


4.95371 


012 










0.65283 


0.28542 


013 










5.67111 


1.87790 


014 










2.06680 


5.46657 


015 










1.78624 


1.14039 


016 












4.75900 


017 












2.68202 


018 












3.51887 


019 












5.54982 


020 












4.35086 


021 












0.478595 



TABLE I: Optimal nearest neighbour concurrences found in FCS, and the corresponding parameters 



V. PROPERTIES OF THE TRANSLATIONAL 
CHAIN MAXIMIZING NEAREST 
ENTANGLEMENT 



Let us investigate the question which properties an in- 
finity translational chain has which maximizes nearest 
neighbour entanglement. Is the obtained chain like an 
ordinary bicycle chain (or Markov chain) , whose links are 
only connected to two neighbouring links but not to the 
next and next-next neighbouring sites as our intuition 
may suggest? 

Lets consider the state of three qubits in a line 



(33) 



and tracing over qubit 2 gives the density matrix for next 



nearest neighbours 

= Y \siS3){tih\ ^Tr{vl^v\vl^pB vt.vivt.,) 

+Tr{vl^v\vl^PB vt,V2VtS^{M) 

If vi is nilpotent, this means that the first term of the the 
last equation vanishes (except of the last one but this has 
just to do with normalization). Thus it reduces the space 
dramatically in which we can vary. We checked that this 
next nearest neighbour entanglement for the parameters 
maximizing the nearest entanglement and found that it 
is zero for all auxiliary dimensions h. Thus we conclude 
that such a chain is an "ordinary bicycle chain" : only two 
neighbouring sites are linked via entanglement. Thisalso 
supports our assumption for a nilpotent generator be- 
cause entanglement concentrates. 

In Fig. 0] we plotted the concurrence and assisted concur- 
rence for the dimension 5 = 2 for different parameters. 
One sees that we can generally find parameters for which 
next nearest entanglement is nonzero, we checked the 
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FIG. 4: (Color online) (a) Here we plot C{pis) for 
ai = 0.6 (blue), 0.42 (red, dashed), 0.3 (green, dotted); the 
line is the maximizing 0i = 0.571859 and the abscissa 
is 01 G [— -I,-!]. Note that we plotted Ai — A2 while 
the concurrence is the maximum of {|Ai — A2|,0}. (b) 
Here the concurrence of assistance Cass ~ pia for ai — 
0.6 (blue), 0.42 (red, dashed), 0.3 (green, dotted) is plotted; 
abscissa G [— ^, ^]. One notices when concurrence in- 
creases, concurrence of assistance decreases and vice versa. 



maximum possible concurrence available for our choice 
of generators, which is below the one of nearest neigh- 
bour entanglement (C13 = 0.169470 for ai = 0.88563 
and 01 — 0.25066 and a nearest neighbour concurrence 
of C12 = 0.270660). When concurrence of next near- 
est neighbour entanglement increases then concurrence 
of assistance of next nearest neighbour entanglement de- 
creases, as it is plotted in Fig^(b). 
As the state of the whole chain is pure, we expect that 
purity has to increase considering the states of more and 
more qubits of the chain. We checked that for all con- 
sidered cases, for example for 6 = 6 Trp\2 < Trpfjs 
is 0.452911 < 0.461722 and for 6 = 9 it changes to 
0.447191 < 0.455342, i.e. far away from 1. Increasing 
dimensionality b increases entanglement and also concen- 
trates the correlation. 

In Table HTl the relevant quantities are summarized. No- 
tice that all quantities behave monotonically with in- 
creasing dimensionality b, except for the length of the 
Bloch vector. The fact that concurrence of assistance is 
bigger than concurrence suggest the presence of multi- 
partite entanglement in the system. 



VI. SUMMARY AND CONCLUSION 



We have addressed the optimization of entanglement be- 
tween nearest neighbours in an infinite translational in- 
variant chain in an increasing set of states. Using finitely 
correlated states with a recursive structure and certain 
reasonable assumptions we have shown that the nearest 
neighbour entanglement almost reaches its conjectured 
upper bound Cw — 0.434467. Our approach allows for an 
explicit calculation of the elements of the nearest neigh- 
bour density matrix maximizing entanglement. We show 
that in a concurrence versus purity (measured by Trp^) 
diagram the obtained nearest neighbour state seems to 
approach the maximally entangled mixed states (MEMS) 
which bound the realizable bipartite states in this dia- 
gram. 

The approach we have adopted has the same roots as 
the DMRG methods. However, instead of using a matrix 
product form for a state of a finite set of qubits the for- 
malism used provides an exact description of the infinite 
chain. The accuracy of the approximation increases with 
the dimensionality b of the auxiliary Hilbert space. 

For dimension 6 = 2 and 6 = 3 we have given a detailed 
analytical treatment of the problem while in higher di- 
mensions we rely on numerical calculations. Apart from 
the investigation of nearest neighbour entanglement we 
have evaluated other properties of the chain. These re- 
sults support the qualitative expectations that due to 
the monogamy of entanglement the increase in the near- 
est neighbour entanglement leads to the decrease in the 
longer distance quantum correlation. Not only concur- 
rence but also the whole nearest neighbour density ma- 
trix tends to reach a given fixed value. This is, however, 
not the case with the state of the whole system. The 
difference of the concurrence and the concurrence of as- 
sistance suggest that some kind of multipartite entangle- 
ment is also present. 

The variational technique utilized in our work may be 
applicable in other similar physical problems. 



Acknowledgments 



The authors wants to thank the "non-local" seminar 
Vienna-Bratislava. B.C. Hiesmayr wants to further ac- 
knowledge the EU-project EURIDICE HPRN-CT-2002- 
00311. M.Koniorczyk acknowledges the support of the 
Hungarian National Grant Agency (OTKA) under the 
contracts Nos. T043287 and T049234, the hospitality of 
Prof. V. Buzek in Bratislava during the first period of 
this project, and the support of the Marie Curie RTN 
network CONQUEST. 



Appendix 



12 



Dimension b 


2 


3 


4 


5 


6 


9 


Cmax(/9l2) 


0.414214 


0.41825 


0.432000 


0.432471 


0.433791 


0.434095 


Cass(pi2) 


0.585787 


0.587251 


0.600000 


0.600131 


0.601204 


0.601442 


A 


0.292893 


0.293626 


0.300000 


0.300066 


0.300602 


0.300721 


a 


0.207107 


0.209126 


0.216000 


0.216236 


0.216895 


0.217048 


c 


0.174155 


0.164125 


0.097378 


0.0925458 


-0.069033 


-0.0575684 


Purity of Trp?2 


0.550252 


0.538009 


0.471242 


0.467748 


0.452911 


0.447191 


Purity of Trpi 


0.646446 


0.639055 


0.598965 


0.597077 


0.58905 


0.586052 


-(Bloch)|2 bTrp|j-l 

1 — b-i 


0.5 


0.607465 


0.536326 


0.554368 


0.519502 


0.521177 



TABLE II: Properties for different dimension b 
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